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1. Introduction 

The meander problem is a simply stated combinatorial question: count the number 
of configurations of a closed non-self-intersecting road crossing an infinite river through 
a given number of bridges. Despite its apparent simplicity, this problem still awaits a 
solution, if only for asymptotics when the number of bridges is large. The problem emerged 
in various contexts ranging from mathematics to computer science In particular, Arnold 
re-actualized it in connection with Hilbert's 16th problem, namely the enumeration of ovals 
of planar algebraic curves , and it also appears in the classification of 3- manifolds . 

Remarkably, the meander problem can be rephrased in the physical language of critical 
phenomena, through its equivalence with a particular problem of Self- Avoiding Walks: the 
counting of the compact foldings of a linear chain. 

Several techniques have been applied to this problem: direct combinatorial approaches 
1^ 0, random matrix model techniques an algebraic approach using the 
Temperley-Lieb algebra and Restricted Solid-On-Solid models 0. Several exact results 
have been obtain on the way for meander-related issues, including exact sum rules for 
meandric numbers 0], the solution of the somewhat simpler irreducible meander problem 
, and the calculation of a meander-related determinant . 

The present paper is dedicated to a more direct enumerative approach and a thorough 
analysis of its results in the spirit of critical phenomena. The meander problem is gener- 
alized to include the case of several non-intersecting but possibly interlocking roads with 
a weight q per road. The corresponding generating functions are analyzed as functions of 
q. In particular, we derive their large q asymptotic expansion in powers of 1/q. 

The paper is organized as follows. In sect. 2, we give the basic definitions of meanders 
and semi-meanders (which correspond to the same problem with a semi-infinite river with 
a source, around which the roads are free to wind), as well as some associated observables 
such as the winding. We further give exact solutions to the meander and semi-meander 
problems at two particular values of q: q = 1, where they reduce to a random walk 
problem, and q = oo, dominated by simple configurations. In sect. 3, we explain how 
to enumerate the semi- meanders for arbitrary number n of bridges, using a fundamental 
recursive construction. After implementation on a computer, this procedure allowed us 
to find the semi-meander numbers with up to ti = 29 bridges. These data are presented 
and then analyzed by a direct large n extrapolation. On the way we also confirm the 
scaling hypotheses borrowed from the theory of critical phenomena. Evidence is found 
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for a phase transition for semi-meanders at a value of q = Qc — 2 between a low-q and a 
large-g regimes, discriminated by the relevance of winding around the source. In sect. 4, we 
show how to use the above data to generate a large-g expansion for most of the interesting 
quantities. This expansion provides an accurate description of the whole q > Qc phase. In 
sect. 5, we analyze the break-down of this expansion, which gives rise to the q < qc phase. 
Sect. 6 briefly describes the small-g expansion of the problem. We gather our conclusions 
in sect. 7. The more technical details are left in appendices. 

2. The meander problem 

2.1. Definitions, observables 

A meander of order n is a planar configuration of a non-self-intersecting loop (road) 
crossing a line (river), through a given number 2n of points (bridges). We consider as 
equivalent any two configurations which may be continuously deformed into each other, 
keeping the river fixed (this is therefore a topological equivalence). The number of in- 
equivalent meanders of order n is denoted by M„. For instance, we have Mi = 1, M2 = 2, 
Ms = 8... More numbers can be found in [ p!0| . 

We stumbled on the meander problem by trying to enumerate the distinct compact 
folding configurations of a closed polymer, i.e. the different ways of folding a closed chain 
of 2n identical constituents onto itself. The best image of such a closed polymer is that 
of a closed strip of 2n identical stamps, attached by their edges, serving as hinges in the 
folding process: a compactly folded configuration of the strip is simply a folded state in 
which all the stamps are piled up on top of one of them. 

Such a compactly folded configuration is easily identified with a meander configuration 
as depicted in Fig.l. Draw a closed line (road) passing though the centers (bridges) of all 
the piled-up monomers, then open one hinge of the polymer (we choose to always open 
the bottom right one) and pull the stamps apart so as to form a straight line: the latter 
is identified with the river, whereas the distorted line becomes the road of the resulting 
meander. 




Fig. 2: The 4 inequivalent foldings of a strip of 3 stamps. The fixed stamp 
is indicated by the empty circle: it is attached to a support (shaded area). 
The other circles correspond to the edges of the stamps. 
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road 

(a) (b) 

Fig. 1: The mapping between compactly folded closed strip of stamps and 
meanders. We display a compact folding configuration (a) of a closed strip 
with 2n = 6 stamps. To transform it into a meander, first draw a (dotted) 
line through the centers of the stamps and close it to the left of the picture. 
Then cut the bottom right hinge (empty circle) and pull its ends apart as 
indicated by the arrows, so as to form a straight line (b): the straight line 
forms the river, and the dashed line the road of the resulting meander. 

When the strip of stamps is open (see Fig.2), we decide to attach the first stamp to a 
support, preventing the strip from winding around it, while the last stamp has a free ex- 
tremal edge. In this case, a slightly generalized transformation maps any compactly folded 
open configuration of (n — 1) stamps to what we will call a semz-meander configuration of 
order n, in the following manner. 




(a) (b) (c) 

Fig. 3: The mapping of a compactly folded configuration of 4 stamps onto 
a semi-meander of order 5. (a) draw a (dashed) curve through the pile of 
stamps and the (shaded) support, (b) pull the free edge of the last stamp to 
form a half-line (the river with a source), (c) the result is a semi-meander 
configuration of order 5, namely that of a road, crossing a semi-infinite river 
through 5 bridges (the source of the river, around which the road is free to 
wind, is indicated by an asterisk). 

As shown in Fig. 3, draw a curve (road) though the (n—l) centers (bridges) of all the 
piled-up stamps, then close this curve across the support (this last intersection is the n-th 
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bridge) , and pull the free edge of the last stamp in order to form a straight half-line (river 
with a source). The resulting picture is a configuration of a road (the curve) crossing a 
semi-infinite river (stamps and support) through n bridges: this is called a semi-meander 
configuration of order n. Note that the road in a semi-meander may wind freely around 
the source of the river, and that consequently the number of bridges may be indifferently 
even or odd, as opposed to meanders. The number of distinct semi-meanders of order n is 
denoted by M^. For instance, we have Mi = 1, M2 = 1, Ms = 2, M4 = 4... More numbers 
can be found in Q and in appendix A. 

Through its compact folding formulation, the semi-meander problem is a particular 
reduction of the two-dimensional self-avoiding walk problem, in which only topological 
constraints are retained. It is therefore natural to define, by analogy with self-avoiding 
walks the connectivity R per stamp and the configuration exponent 7 which determine the 
large n behavior of the semi-meander numbers as follow^ 

Mn ~ c — (2.1) 

The connectivity R may be interpreted as the average number of possibilities of adding 
one stamp to the folded configurations. The exponent 7 is characteristic of the (open) 
boundary condition on the strip of stamps. 




(a) (b) 

Fig. 4: The "end-to-end distance" of the folded strip of stamps (a) is the 
number {w = 1 here) of stamps to be added to the strip (the added stamp 
is represented in dashed line), so that the new free end (empty circle) is in 
contact with the infinity to the right. This coincides with the "winding" of 
the corresponding semi-meander (b), namely the number of bridges to be 
added if we continue the river to the right of its source (dashed line). 



That the semi-meander numbers M„ actually have these leading asymptotics may be proved 
by deriving upper and lower bounds on R. See Q for further details. 
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A natural observable for self-avoiding walks is the end-to-end distance. The corre- 
sponding notion for a compactly folded open strip of stamps is the "distance" between 
the free end of the strip and, say the support. This distance should also indicate how far 
the end of the strip is buried inside the folded configuration. It is defined as the minimal 
length w of a, strip of stamps to be attached to the free end, such that a resulting folding 
with n — 1 + w stamps has its free end outside of the folding, namely can be connected to 
the infinity to the right of the folding by a half-line which does not intersect any stamp. 
Indeed, the infinity to the right can be viewed as the nearest topological neighbor of the 
support, hence w measures a distance from the free end of the strip to the support. This 
is illustrated in Fig. 4(a), with n = 5 and w = 1. In the semi-meander formulation (see 
Fig.4(b)), this distance w is simply the winding of the road around the source of the river, 
namely the number of bridges to be added if we continue the river to the right of its 
source. By analogy with self-avoiding walks, we expect the average winding over all the 
semi-meanders of order n to have the leading behavior 

= aT E ^ ~ (2.2) 



semi — meanders 



where u is some positive (end-to-end) exponent < < 1, as w is always smaller or equal 
to n. 

In this language, a meander of order n is simply a semi-meander of order 2n with 
winding w = 0. By analogy with closed (as compared to open) self-avoiding walks, we 
expect the asymptotics 

~ c ^ (2.3) 

where the connectivity per bridge R is the same as that for semi-meanders ( |2.1D , R = R, 
and the configuration exponent a 7^ 7 is characteristic of the closed boundary condition 
on the strip of stamps. 

In the following, we will mainly focus our study on the semi-meander numbers. 

2.2. Arches and connected components 

Any semi-meander may be viewed as a particular meander by opening the semi-infinite 
river as indicated by the arrows on Fig. 5. In the process, the number of bridges is doubled, 
hence the order is conserved. The resulting meander however is very peculiar. Note that 
in general a meander is made of an upper (resp. lower) configuration consisting of non- 
intersecting arches (arcs of road) connecting the bridges by pairs above (resp. below) the 
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Fig. 5: A semi-meander viewed as a particular meander: the semi-infinite 
river must be opened up as indicated by the arrows. This doubles the number 
of bridges in the resulting meander, hence the order is conserved (ti = 5 
here). By construction, the lower arch configuration of the meander is always 
a rainbow arch configuration of same order. 

river. In the present case the lower configuration is fixed: it is called the rainbow arch 
configuration of order n (the bridge i is connected to the bridge {2n — i + l),z = l,2,...,n). 
On the other hand, the upper arch configuration may take any of the Mn values leading 
to semi- meanders of order n. 
There are however 

n!(n -I- 1)! 

distinct arch configurations of order n [0, as is readily proved by recursion (cn+i = 
YjQ<j<n^3^ri-j, with cq = 1, heuce ci = 1, C2 = 2, C3 = 5, C4 = 14,...: the are 
called the Catalan numbers). Hence not all upper arch configurations, once supplemented 
by a lower rainbow arch configuration of same order, lead to an opened semi-meander 
(M„ < Cji). This is because, in general, the corresponding object will have k > 1 con- 
nected components: we call it a semi-meander of order n with k connected components. 
Indeed, if the river is folded back into a semi-infinite one, we are simply left with a col- 
lection of k possibly interlocking semi-meanders of respective orders rii, 712,..., n^, with 
ni -f 77.2 + ... + nfc = n. We always have 1 < /c < n, and k = n only for the superposition of 
an upper and a lower rainbow configurations, leading to n concentric circles. We denote by 

~ (k) 

Mn the number of inequivalent semi- meanders of order n with k connected components. 
In particular, we have Mn"^ = Mn and Mn"^^ = 1 for all n. 

^ (k) 

The direct numerical study of the asymptotics of the numbers MA turns out to be 
delicate, as the natural scaling variable of the problem is the ratio x = k/n, which depends 
on n and takes only a discrete set of values. To circumvent this problem, we will study 
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the generating function mn{q) for these numbers, also referred to as the semi-meander 
polynomial. 

n 

mniq) = $^?'M^'^ (2.5) 

k=i 

This quantity makes it possibfe to study the large n asymptotics of the Mn''^ in a global way, 
by use of extrapolation techniques for all real values of q. The semi-meander polynomial 
( ^■5| ) may be viewed as the partition function of a statistical assembly of multicomponent 
semi-meanders of given order n, with a fugacity q per connected component. As such, it 
is expected to have an extensive large n behavior, namely 

~ -c{q) ^ (2.6) 

where R{q) is the partition function per bridge, 7(13') is a possibly varying exponent and 
c{q) a function independent of n. For g ^ (/c = 1), we must recover the connected 
semi-meanders, namely that rfi„,{q)/q M^, i.e. 

R{q)^R 7(9) -^7 ciq)/q^c (2.7) 

(c.f. (|2.1| )). The notion of winding is well-defined for multi-component semi-meanders as 
well, as the sum of the individual windings of each connected component, namely the total 
number of times the various roads forming the semi-meander wind around the source of 
the river. Therefore we define 



nT-niq) 

^ ' multicomp. 

semi — meanders 

where u{q) is the generalized winding exponent for multi-component semi-meanders, sat- 
isfying < z^(q') < 1. 

Analogously, we define multi-component meanders of order n, as configurations of k 
non-intersecting roads {1 < k < n) crossing the river through a total of 2n bridges, and 
denote by their number. We also define the meander polynomial 

n 

mniq) = E^' (2.9) 
k=i 

This is nothing but the restriction of (p.5| ) with n 2n, to semi-meanders with zero 
winding w = 0. We therefore expect the asymptotics for large n 

mM - c{q)^l^ (2.10) 



In this estimate, the partition function per bridge R{q) is expected to be identical to that 
of semi-meanders R{q) only if the winding is irrelevant, namely if u{q) is strictly less than 
1 

R{q) = R{q) iff i^{q) < 1 (2.11) 

Otherwise, the fraction of semi-meanders with zero winding may be exponentially small, 
and we only expect that R{q) < R{q) if i^iq) = 1- 

2.3. Exact results for large numbers of connected components (q = oo) 
For very large q, we simply have 



mn{q) 



q 



(2.12) 



as the meander polynomial is dominated by the k = n term, corresponding to the unique 
semi-meander of order n made of n concentric circular roads, each crossing the semi-infinite 
river only once. This semi-meander will appear as the rightmost object in the n-th line of 
the tree of Fig. 8. The winding of this semi-meander is clearly w = n, hence we have, for 
q ^ oo 

R{q)^q 7(9)^0 c{q)^l iy{q) ^ 1 (2.13) 

As to meanders, the only way to build a meander of order n with the maximal number 
n connected components is that each component be a circle, crossing the river exactly twice. 
This is readily done by taking any upper arch configuration and completing it by reflection 
symmetry w.r.t. the river. This leads to Mn""* = Cn (c.f. ( p.4|) ) meanders with n connected 



components. By Stirling's formula, we flnd that when q 
behaves as 

rriniq) ~ q"^ 

1 (2v/?)''^ 



00 the meander polynomial 



(2.14) 



n n 



3/2 



hence, when q 00 



R{q) ^ 2v^ a{q) 3/2 c{q) ^ 1/v^ 



^2.15) 



This conflrms the abovementioned property (|2.11|) that R{q) < R{q) when i'{q) = 1, as 
2^ < Q for large q. 



2.4- Exact results for random walks on a half-line (q = 1) 



When q = 1 in ( p.5| ), rn^(l) simply counts all the multi-component semi-meanders, 
irrespectively of their number of connected components. This simplifies the problem dras- 
tically, as we are simply left with a purely combinatorial problem which can be solved 
exactly. The multicomponent semi-meanders are obtained by superimposing any arch 
configuration of order n with the rainbow of order n, hence 

1 4'^ 

m„(l) = Cn ^^77? (2.16) 



by use of Stirling's formula for large n. This gives the values 

i?(l)=4 7(1) = 3/2 c(l) = l/v^ (2.17) 



The study of the winding at = 1 is more transparent in the formulation of arch 
configurations of order n as random walks of 2n steps on a semi-infinite line. For each arch 
configuration of order n, let us label by 1, 2,...,2n — 1 each segment of river in-between two 
consecutive bridges, and the leftmost semi-infinite portion, 2n the rightmost one. Let 
h{i), i = 0, 1, 2?7. denote the number of arches passing at the vertical of the corresponding 
segment i. By definition, h{0) = h{2n) = 0. More generally, going along the river from 
left to right, we have h{i) = h{i — 1) -|- 1 (resp. h{i) = h{i — 1) — 1) if an arch originates 
from the bridge i (resp. terminates at the bridge i). 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



Fig. 6: A walk diagram of 18 steps, and the corresponding arch configuration 
of order 9. Each dot corresponds to a segment of river. The height on the 
walk diagram is given by the number of arches intersected by the vertical 
dotted line. 
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The function h satisfies h{i) > 0, for all i, and may be interpreted as a "height" 
variable, defined on the segments of river, whose graph is nothing but a walk of 2n steps 
as shown in Fig. 6. This may be seen as the two-dimensional extent of a brownian motion 
of 2n steps on a half-line, originating and terminating at the origin of the line. This 
interpretation makes the leading behavior ~ 2^*^ of ( p.l6|) clear: it corresponds to the 
2 possible directions (up or down) that the motion may take at each step. The exponent 
3/2 in ( 2.161 ) is characteristic of the boundary condition, namely that the motion is closed 



and takes place on a half-line (other boundary conditions would lead to different values of 
7, e.g. for a closed walk on a line, we would have a behavior (^2^^) ~ l?"^ j ^Jn). 

In this picture, the winding is simply given by the height w = h{n) of the middle point. 
Let us evaluate more generally the average height of a point i over the arch configurations 
of order n. It is given by 

(/i(z))„ = —\2hA^,,{h) (2.18) 
c 

where A„ ^(/i) denotes the number of arch configurations of order n such that h{i) = h. A 
simple calculation p shows that 

..(.^((4)-(^VO)(C^-"4)-C-'^VO) 

as the An^i{h) walks are simply obtained by gluing two independent walks of i and 2n ~ i 
steps linking the origin to the height h. 

In the case of the winding, w = h{i = n), (|2.18| ) leads to a more compact formula, 
according to the parity of n 



n = 2p: {w)2p = — 1 

C2p 

r2p\ /2p+i\ 

n = 2p + l: {w)2p+i = 2 ^^^^ ^ ' -1 

C2p+1 



(2.20) 



For large n, this gives the following expansion 



{w)^ = 2./^ - 1 + + 0(l/n3/2) (2.21) 



irrespectively of the parity of n. This implies that 

u(q = 1) = 1/2 (2.22) 
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This is the weU-known result for the Brownian motion, for which the extent of the path 
scales like n^/^ for large n. It is instructive to note that, thanks to ( |2.21| ), the observable 
w + 1 is less sensitive than w to the finite size effects at q = 1. This will be useful in 
the forthcoming numerical estimates for arbitrary q where we observe that the numerical 
extrapolations are improved by considering w + 1 instead of w. Using ( |2.19| ), we may now 
compute the probability distribution Pn{w) for an arch configuration of order n to have 
winding h{n) = w, which takes for large n the scaling form 

Pr^iw) = -A„,nH ~ yV/fy^l (2-23) 



Cn ' {w)n \{w) 

with a scaling function / independent of n for large n, readily obtained by use of Stirling's 
formula, upon writing w = ^^njix ^ for large n. This gives 



/(O = -^ee-^^^' (2.24) 



for all ^ > 0. 

For general position i ^ n, M^e find, by a saddle point evaluation of the sum ( |2.18| ), 
that the average profile of arch configurations is a "Wigner" semi-circle 



{h{i))^ ~ 2^-^x{2-x) (2.25) 

when expressed in the scaled position x = i/n, 0<x<2. 

The meanders of order n are the semi-meanders of order 2n with winding w = h{2n) = 
0. They are therefore built as the juxtaposition of two independent walks of length 2?!. 
Hence 

m,(l) = (cn)' ^ (2.26) 

or, in other words 

R{l)=R{l)=A a(l)=3 c(l) = l/7r (2.27) 

This is again in agreement with ( |2.11|) , as z^(l) = 1/2 < 1, i.e. the winding is irrelevant at 
q = l. 
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3. Exact enumeration and its analysis 

In this section, wc present results of an exact enumeration of Mn'^^ for small n (n < 29), 
and analyze their large n extrapolation. The enumeration is performed by implementing 
on a computer a recursive algorithm which describes all the semi-meanders up to some 
given order. Clearly, the complexity is proportional to the Catalan numbers (c„ ~ 4"^) 
hence the limitation on n. 

3.1. The main recursion relation 

The subsequent numerical study relies on the exploitation of the following recursion 
relation generating all the semi-meanders of order (n -|- 1) from those of order n. 




Fig. 7: The construction of all the semi-meanders of order n+1 with arbi- 
trary number of connected components from those of order n. Process (I): (i) 
pick any exterior arch and cut it (ii) pull its edges around the semi-meander 
and paste them below. The lower part becomes the rainbow configuration 
TZn+i of order n+1. This process preserves the number of connected com- 
ponents k ^ k. Process (II): draw a circle around the semi meander of order 
n. This process adds one connected component k ^ k+1. 

We start from any semi-meander of order n with k connected components, in the 
open- river picture. We may construct a semi- meander of order {n-\-l) in either following 
way (denoted (I) or (II)), as illustrated in Fig. 7 

(I) Pick any exterior arch, i.e. any arch with no other arch passing above it. Cut 
it and pull its ends all the way around the others (in order to add two bridges), and 
reconnect them below, by creating an extra concentric lower arch for the rainbow. In this 
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process, we have n — > n + 1, but the number of connected components has not changed: 
k ^ k. Another way of picturing this transformation is the following: one simply has 
pulled the exterior arch all the way around the semi-meander and brought it below the 
figure, creating two new bridges along the way. As no cutting nor pasting is involved, the 
number of connected components is clearly preserved. 

(II) Draw a circle around the semi-meander. This adds a lower concentric semi-circle 
which increases the order of the rainbow to (n+l), and also adds one connected component 
to the initial semi- meander k ^ k + 1. 

These two possibilities exhaust all the semi- meanders of order {n+ 1), as the trans- 
formation is clearly invertible, by pulling back up the lower external arch of the rainbow. 
Note that by construction, there are as many possibilities for the process (I) as exterior 
arches, and the transformation is therefore one-to-many. 




Fig. 8: The tree of semi-meanders down to order n = 4. This tree is 
constructed by repeated applications of the processes (I) and (II) on the semi- 
meander of order 1 (root). We have indicated by small vertical arrows the 
multiple choices for the process (I), each of which is indexed by its number. 
The number of connected components of a given semi-meander is equal to 
the number of processes (II) in the path going from the root to it, plus one 
(that of the root). 

We may now construct a tree of all the semi-meanders, generated recursively from 
that of order 1 (root), as displayed in Fig. 8. Note that we have adopted the open-river 
formulation to represent them. 

Keeping track of the connected components, this translates into the following relation 
between the semi-meander polynomials 

fhn+i{q) = fhn{q) {ext. arch.) n{q) + qfhn{q) (3.1) 
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where we denoted by (ext. arch. the average number of exterior arches in a semi- 
meander of order n, weighed hy q'^, k its number of connected components. In ( |3.1| ), the 
first term corresponds to all the processes (I), whereas the second term corresponds to (II). 
Taking the large n limit in ( |3.1| ), this permits to interpret 



R{q) - Q = (ext.arch.)oo(?) (3.2) 

as the limit when n — > oo of the average number of exterior arches in semi-meanders of 
order n, weighed by an activity q per connected component. 
In particular, when g oo, we have 

Riq) -q-^1 (3.3) 

as the limiting semi-meander, made of concentric circles, has only one exterior arch. This 
is a refinement of the large q estimate in ( p.l3|) . When g = 1, we find an average of 



(ext.arch.)oo(l) = R{1) - I = 3 (3.4) 

exterior arches in arbitrary arch configurations of large order 0. When q = 0, the con- 
nected semi- meanders, with k = 1, are obtained through repeated action of the process (I) 
only. This restricts accordingly the tree of Fig. 8. In that case, the partition function per 
bridge 

R = R(0) = (ext.arch.)oo(O) (3.5) 

coincides with the average number of exterior arches in connected semi-meanders, for large 
n. 



3.2. Exact enumeration and large n extrapolation 

In appendix A, we give an archetypical example of the programs we have implemented 
to compute the semi-meander numbers and various observables. We have computed the 
numbers MA for 1 < k < n < 27, the numbers MA for n < 29. To investigate the 
winding of semi- meanders, we have also computed the numbers Mn'^'' (w) of semi-meanders 
of order n, with k connected components and fixed winding w for 1 < k < n < 24 (encoding 
in particular the meander numbers Mn'^'' = M2^{0)), and the semi- meander "profile" 

E (3.6) 

semi — meanders 
order n ,k c.c. 
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for all positions < z < 2n. Some of these numbers can be found in appendix A. 

After gathering these numbers into generating functions of q, it is possible to perform 
large n extrapolations at fixed q, for the quantities R{q), R{q), 7(?), c{q), c{q) and 

u{q). 

The general extrapolation scheme is the following. Suppose an observed quantity 
has the following large n expansion 

X. = E^ + 0(l/n-+^) (3.7) 

A;=0 

Then we get a best estimate of the large n limit xq by iterating p times the difference 
process (A/)(n) = f{n + 1) — f{n) on the function f{n) = n^Xn, with the result 

— n^X^ = xo + 0{l/nP+^) (3.8) 

This gives perfect results for the Catalan numbers (i.e., q = I) using X^ = Log(c„_|_i/c„). 
This turns out to extend to a whole range of g's in a neighborhood of 1. For instance, 
Log R{q) is extrapolated using X^ = Log 1^771^+1 (<?) / ^n-i (?) • 

The results for R{q) and R{q) are displayed in Fig. 9. The two functions are found to 
coincide in the range < q < qc with qc — 2, and to split into R{q) > R{q) for q > qc- 
As explained before, the comparison between R{q) and R{q) determines directly whether 
u{q) is 1 or not. The result of Fig. 9 is therefore the signal of a phase transition at q = qc 
between a low-q regime where the winding is essentially irrelevant {iy{q) < 1) and a large-g 
phase with relevant winding {iy{q) = 1). 

This is compatible with the direct extrapolation for u{q) displayed in Fig. 10, which 
is however less reliable in the region around q = 2, due to its sub- leading (and probably 
discontinuous) character. 

The configuration exponent for semi- meanders 7(5) is represented in Fig. 11, for two 
different orders in the extrapolation scheme ( p.8| ). The extrapolation proves to be stable 
for < < 2. For g > 2, it develops oscillations around a mean value, estimated to vanish 
(^(qi) 0) for q large enough. For simplicity, we chose not to represent the functions a{q), 
c{q), c{q). The coefficient c{q) develops a discontinuity at the transition q = 2. On the 
other hand, the functions pertaining to meanders only {a{q) and c{q)) do not display any 
transition aX q = 2. 
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Fig. 9: The functions R{q) and R{q) for < < 4 as results of large n 
extrapolations. The two curves coincide for < < 2 and split for q > 2 
with R{q) > R{q)- Apart from the exact value -R(l) = -R(l) = 4, we find the 
estimates R{0) = 3.50(1), R{2) = 4.44(1), R{3) = 4.93(1) and i?(4) = 5.65(1). 

3. 3. Scaling functions 

By analogy with critical phenomena, in addition to the scaling behaviors ( p.6|) , ( |2.1C| ) 
and (|2.8| ) involving the critical exponents 7(1?), a{q) and i^{q), we expect to find more 
refined scaling laws involving scaling functions. A particular example of such scaling 
functions has been derived for g = 1 ( p.23| ), for the probability distribution Pn{w) of the 
winding w among arch configurations of order n. It involves the scaling function ( |2.24| ). 
For q = we expect the same behavior for the corresponding probability distribution 



of winding w among connected semi-meanders of order n. We expect the scaling behavior 




(3.9) 




(3.10) 
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Fig. 10: The winding exponent for < g < 8, as obtained from a 
large n extrapolation. We observe a drastic change of behavior between low 
g's and large g's, with an intermediate regime where the extrapolation fails, 
hence is not reliable. The dashed line indicates a possible scenario for the 
exact function iy{q), compatible with a transition at Qc — 2. Apart from the 
exact value = 1/2, we read i/(0) = 0.52(1). 

This is precisely what we observe in Fig. 12, where we plot {w + l)n(O) Pn^\w) as a 
function of the reduced variable ^ = {w + 1)/ {w + l)n(O) for different values of n. Indeed, 
as already explained in the q = 1 case, we have taken the variable {w + 1) instead of w to 
improve the convergence. All the data accumulate on a smooth curve, which represents the 
scaling function f^'^\^). The shape of this function is reminiscent of that of the end-to-end 
distribution for polymers. By analogy, we expect a certain power law behavior for small ^ 



where 6 satisfies the relation 



obtained by identifying 



to 



- e (3.11) 



a -7 = iy{l + e) (3.12) 



^in^(O) oc ^f^'^ (^) (3.13) 



^ = ^ocn- (3.14) 
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Fig. 11: The configuration exponent 7(g) for < g < 4, from two different 
large n extrapolations. Apart from the exact value 7(1) = 3/2, we estimate 
7(0) ~ 2. 
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Fig. 12: Plot of {w + 1)^(0) Pi^\w) function of the reduced variable 
^ = (w+l)/ {w+l)n{0) for n — 2,3, 24. The points accumulate to a smooth 
scaling function f^^\^). The erratic points correspond to small values of n, 
which have not reached the asymptotic regime. 
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For large ^, we expect a behavior f^^\^) ~ exp(— const. with a possible Fisher-law 
behavior 6 = 1/(1 — u). The observed function of Fig. 12 is compatible with these limiting 
behaviors, although we cannot extract reliable estimates of the exponents 9 and 5. 

As we already did in the case oi q — 1 ( |2.25| ), we can study the average profile of 
semi- meanders 

m)niq) 



—pix = i/n]q) 
n 



(3.15) 



involving a scaling function q) of the variable x, with < x < 2 for each value of q 
(with the appropriate normalization such that j p = 1). For instance, we have seen in 
eq.(|3|) that p(x; 1) = {2/t:)^x{2 - x). 
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Fig. 13: Semi-meander average profiles for g = 0, 1, 2, 4, and 1 < n < 24, 
as functions of the reduced variable x = i/n. For g = 1, we also represented 
the exact large n Wigner semi-circular limit p{x; 1). For g = 4, we also 
represented the large n and q piecewise-linear limit p{x;oo). 



We have represented in Fig. 13 these profiles for several values of q. Again, the points 
accumulate on smooth limiting curves p{x;q). We observe a first change of behavior at 
q = 1 between a q < 1 regime with a negative cusp at x = 1 and a q > 1 regime with a 
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positive cusp, separated by the Wigner semi-circle, with no cusp at g = 1. For large g, 
p{x; q) tends to the limit p{x; oo) = 1 — |1 — x| corresponding to the unique semi-meander 
made of n concentric circles, which satisfies h{i) = i for < z < n and h{i) = 2n — i for 
n < i < 2n. For small x, we expect a power law behavior of the form 

p{x;q) ~ (3.16) 

where we identify the exponent f>{q) = i^{q) from h{l) = 1 and the fact that h{j) ~ 



4. Large q asymptotic expansions 

In the previous section, we have observed two regimes for the semi-meander polyno- 
mials, namely a \ow-q regime in which the winding is irrelevant and a large-g regime where 
the winding is relevant, separated by a transition at a value of q = qc — 2. On the other 
hand, we have already exhibited an exact solution of the problem at g = oo (|2.13| ), and a 
first correction thereof for large q in ( |3.3| ). It is therefore tempting to analyze the large q 
phase by a systematic expansion in 1/q. This is performed in the following section, where 
R{q) is expanded up to order 19 in 1/q, and 7(5) is found to vanish identically throughout 
the large q regime. In the subsequent section, we compute the large q expansion of the 
average winding in semi-meanders, and we find u{q) = 1 identically in this regime. 

4.1. Large q asymptotic expansion of the semi-meander polynomial 

In this section we derive the large q expansion of the semi-meander polynomial rhn{q) 
of eq.(|2.5D as 



rnniq) = (M^^ + — + + " " O (4-1) 

q q^ 

involving the semi- meander numbers in the form MA , /c = 0, 1, 2,... Remarkably, these 
numbers, for arbitrary n > 2/c — 1, are polynomials of n, which furthermore exhibit some 
special structure allowing for an explicit large q expansion of R{q). 

The section is organized as follows. We first derive the polynomial form of the MA , 
valid for n > 2A; — 1, together with the corrections to be added for smaller n's. 

The re-exponentiation of rhn{q) in the form ( p.6|) induces strong constraints on the 
polynomials MA , which allow for their complete determination up to /c = 18, out of 
their first values for small n, which were enumerated exactly up to n = 27. 
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Fig. 14: Semi-meanders with many connected components, {a) k = n con- 
nected components; there are n circles. (6) k = n — 1 connected components; 
there are (n — 2) circles and one "kidney" . [c) k = n — 2 connected compo- 
nents; there are respectively (c)i two disjoint kidneys and (n — 4) circles; (0)2 
two kidneys included in one another and (n — 4) circles; (0)3 and (0)4 one 
"spiral" and (n — 3) circles. 

For starters, let us first compute the numbers Mn"^ for /c = 0, 1, 2. 

As we already mentioned, the leading term Mn"^^ = 1 in the expansion counts the 
unique semi-meander of order n with n connected components made of concentric circles 
only, and which we refer to as the leading semi-meander (c.f. Fig. 14 (a)). This yields the 
first polynomial 

Po{n) = M^) = 1 (4.2) 

for all n > 0. 

The sub-leading term is made of the 

pi(n) = M^"-^) = n-1 (4.3) 

"kidney"-type perturbations of the leading semi-meander, displayed in Fig. 14 (b). 

The next-to-leading term consists of the semi-meanders of order n with (n — 2) con- 
nected components, which may be obtained as follows. 

(i) A first possibility consists in taking two "kidney" -type perturbations of the leading 
semi-meander (see Figs. 14 (c)i and (0)2), which are either 

(il) (Fig. 14 (c)i) disjoint, hence a total of (n — 2){n — 3)/2 choices for the positions 

of two kidneys. 
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(i2) (Fig. 14 (0)2) included in one another, hence a total of {n — 3) choices for the 
position of the double kidney, or choice if 71 = 2. 

(ii) The second possibility is a larger "spiral" -type perturbation of the leading semi- 
meander, with a total of {n — 2) available positions, and there are two such perturbations 
(see Figs. 14 (0)3 and (0)4). 

Summing up all these contributions givesi 



hence a polynomial 



/ N n'^ + n- 8 , 
P2{n) = (4.5) 



We see here the first appearance of a correction for small n, in that an extra boundary 
term (5n,2) has to be added at n = 2, to recover the fact that M^^^ = 0. 

More generally, let us consider the number Mn"^ for large n and finite k. The 
corresponding semi- meanders are obtained in the tree of Fig. 8 by applying to the root k 
times the process (I) and (n — 1 — /c) times the process (II). This gives ("^^) ~ /k\ 
possible choices {k « n). These choices however are not completely independent. 

Recall that the process (I) may be applied to any exterior arch of the semi- meander. 
In the situation where k << n, the semi-meander will most probably have only one exterior 
arch (generated by the last process, most probably of the type (II)), and there will be only 
one choice for (I), creating a kidney. If two or more processes (I) are applied consecutively, 
the number of exterior arches may grow, and yield more semi-meanders (e.g. the appli- 
cation of two consecutive processes (I) yields three possibilities: two included kidneys, or 
any of the two spirals). Such an effect is however sub- leading, as the number of choices of 
two or more consecutive processes (I) grows at most like 'n^~^ . Collecting all these extra 
combinatorial factors permits, like in ([4.4|), to write the number = pk{n) as a 



polynomial of degree k in n for large enough n, and fixed k, with moreover a leading term 
n'^/kl. 

For smaller n, the above consecutive choices are affected by the boundaries of the 
tree. Recall for instance that the double included kidneys of Fig. 14 (0)2, obtained by two 
successive applications of (I) on two concentric circles, may only exist for n > 4. We 



Throughout this section, we use the fact that M^°^ = 5n,o, i-e. there are no semi-meanders 
with zero connected component, except for the vacuous semi-meander of order 0. 
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have therefore needed to specify that their number (genericaUy equal to (n — 3)) vanishes 
for n = 3 (granted) and n = 2, the latter resulting in a boundary correction 5^,2- More 
generally, expressing that the combinatorial expressions found are only valid for large 
enough n's will translate into boundary terms, studied in detail in appendix B. We simply 
quote the result here, valid for all n> k 

k-2 

M^n-^^ = Pkin) + J2 fif^ K2k-2-j (4.6) 

where the fl^j'^^ are some positive integers. In particular, ( [4.6| ) shows that the formula 
Mn"^ is a pure polynomial Pkin) with no corrections as soon as n > 2/c — 1. This 

(k) 

property is derived in appendix B, where the first correction Hq = Ck-i is also obtained. 

(k) 

The precise determination of the p^'s and the could be in principle achieved 
directly by pursuing the above method used for /c = 0, 1, 2. However the complexity of this 
program is comparable to that of the exact enumeration of the semi-meander numbers. 
Instead we can guess the coefficients of Pk{n) by matching our data for Mn"^ with the 



form ( [4.6| ). This can be pushed further by exploiting the re-exponentiation property of 
rhn{q), which implies relations between the coefficients of the Pk{n), as discussed now. 

The property ( [4.6|) must be reconciled with the large n behavior of rfin{q) (|2.6| ), namely 
that for n sufficiently large 

Logfhniq) = nLog R{q) - 'J {q)Log n + Log c{q) + o{l) (4.7) 

Such an expansion is valid for all q, and we can in particular study it for large q. On the 
other hand, up to any order k in the 1/q expansion, and by choosing n > 2/c — 1, we may 
also write 

Log^ = Log(l + ^ + ^ + ... + ^ + 0(^)) 
q-^ ^ q q^ q^ \k+i^J 

= E^EH~ E PkAn)PkAn)...Pk,{n) + Oi-j^) 

m=l ^ j = l fei + ...+fej=m y 

fej>l 

In this expansion, the coefficient of 1/q'^ is a polynomial of n, as a sum of products of 
polynomials of n. Comparing this with an expansion of (|4.7| ) in 1/q, we see that its degree 
is at most 1. Therefore, there exist two sequences of coefficients {ak,Pk), such that 

Pkin) = akn + /3k + — : — ^ Pfci (?^)Pfc2 (^)---Pfc,- (n) (4.9) 

j=2 ^ ki + ...+kj=k 
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with the correspondence 



Log R{q) 
Log c{q) 



k>i ^ 



qk 

k>i ^ 



(4.10) 



Moreover, there can be no Logn term in the expansion of (|4.7|) , hence the remarkable 
result 

7(9) = (4.11) 
This result is expected to hold as long as the corrections to the polynomial behavior of the 



(n-fc) 



are negligible. As we will see, this condition defines precisely the large q phase 



q > qc- Therefore the exponent 7(5) vanishes identically over the whole phase q > qc- In 
view of, say, the exact value 7((j = 1) = 3/2, this property cannot persist in the small 
q < qc phase. This is not surprising since the first correction fx^^^ = Ck-i (for k = n/2 + 1) 
implies an additional power law correction of the form 1/n"^/^. 

For n > 2k — 1, (|4.9| ) is a quasi-recursion relation for the polynomials pk, hence for the 
semi- meander numbers M^T ■ This relation is exploited in appendix C to generate from 
our numerical data the polynomials Pk{x) for < /c < 18, together with their corrections. 
Using these polynomials, it is now straightforward to read the functions R{q) and c{q) 
using ( |4.10| ), with the result 
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(4.12) 



It is interesting to compare the result of these large q expansions to the previous direct 
large n extrapolations of sect. 3. As far as R{q) is concerned, we find a perfect agreement 
for the values g > 2, down to = 2, where we find -R(2) ~ 4.442(1) using ( [4.12| ), in perfect 
agreement with the previous estimate. The precision of ( f4.12| ) increases with q, leading to 
far better estimates than before: ^(3) ~ 4.92908(1), ^(4) ~ 5.6495213(1)... 
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As to 7(0'), our prediction that 7(g) = for all > 2 is compatible with the previous 
extrapolation of Fig.ll, where this value is represented in dashed line. We therefore expect 
7((j) to have a discontinuity at g = 2, where it goes from a non-zero 7((j = 2~) value to 
zero. This will be confirmed by the forthcoming analysis of the low-g phase in sect. 5. 

Large q asymptotic expansion of the semi-meander winding 

In this section, we first examine the contribution of the circles to the winding of semi- 
meanders in the large q phase. It turns out to be of the order n throughout this phase, 
implying that the winding exponent u{q) = 1 for all q > qc- In a second step, we compute 
the large q expansion of the average winding by techniques similar to those of the previous 
section, showing that the circles contribute only for a finite fraction of the winding in the 
large q phase. This study will single out the value qc = "2- with a very good precision. 

To enumerate the total number of circles in order n semi-meanders, we simply have 
to count the semi-meanders with a marked circle, in one-to-one correspondence with pairs 
of semi-meanders of total order (n — 1) (since the marked circle separates the original 
meander into two disconnected pieces, its inside and outside). Hence, the average number 
of circles in semi-meanders is given by 

{circ.)n{q) = q— (4.13) 

valid for all q and n, with the convention that rfio{q) = 1. Using the large n asymptotics 
^n{,<l) ~ c{q)R{q)'^, with 7(5) = throughout the large q regime, we find 

{circ.)M -nq^ (4.14) 
R[q) 

The average number of circles therefore grows like n which in turn implies that 

v{q) = 1 q > qc (4.15) 

throughout the large q regime, since each circle contributes 1 to the winding and clearly 
{w) > (circ). Note that the above argument relies crucially on the fact that 7(5) = 0, 
and thus cannot be applied to the small q regime. Indeed, for = 1, ( [4.13| ) simply gives 
(circ.)^(q' = 1) = 1 for all n, hence a very different behavior. 
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The quantity a{q) = lim (circ.)j^(g)/n is simply given by the Taylor expansion of 

n — >oo 

qc{q)/R{q) 
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(4.16) 



Let us now turn to the large q expansion of the average winding of semi-meanders 
( p.2| ). This requires a refined study of the semi-meander numbers '^\'>^) with fixed 

winding w, which display a similar polynomial structure as the The study of the 



corresponding generating function is presented in appendix D and leads to 



{w)niq) = Xiq)n + n{q) 



(4.17) 



where the coefficients X{q) and fj,{q) have the following large q expansions up to order 14 
in 1/q 
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(4.18) 



It can be checked directly that X{q) > a{q) hence the circles only contribute for a 
finite fraction of the total winding. The plots of the functions X{q) and a{q) are displayed 
in Fig. 15. Remarkably, both coefficients seem to vanish at the same point q = 2 with an 
excellent precision. Since these coefficients must be positive, we deduce that our large q 
formulas break down for q < 2. We interpret this as yet another evidence of the drastic 
change of behavior of the average winding (w), which is no longer linear in n below qc, and 
we find qc = 2 with an excellent precision. 

The actual break-down of the large q phase is studied more systematically in next 
section. 
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Fig. 15: The series X{q) ( [4.18|) and a{q) ([4 .161 ) of 1/q up to order 14 and 18 



respectively, for 1 < q < 8. Both curves seem to vanish at q — 2. 
5. The break-down of the large q expansion for q < qc 

The properties of the large q phase are intimately based on the polynomial structure of 
the numbers MX ■ In this section, we explain the break-down of this phase by the precise 



structure of the non-polynomial corrections ( [4 .61 ) to this behavior. The phase transition 
occurs when these corrections become dominant. The detailed analysis of these corrections 
shows a strong resemblance between the low q phase and a meander (zero- winding) regime. 
The corrections to the polynomial behavior of the MA are gathered in the function 

n n 

= 5^(M(-'=)-p.(n))g--^ = Yl ^2t-2-n?^"' (5-1) 

k=0 fc=[(n+2)/2] 

The numbers Ji2k_2_n listed for n/2 < /c < n = 2, 3, 9 in appendix C. With the only 
difference that we are now dealing with rational fractions of n instead of polynomials, we 
can carry the same large n expansion as in sect. 4.1, to extract the large n asymptotics of 
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UniQ), with the result, according to the parity of n 



{2n - 1)^1(9) 



(5.2) 



Using the results of appendix C, we get the following large q expansions (incorporating 
the large n asymptotics of Cn ~ 4'^ /{y/nn^^'^)) 
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(5.3) 



Due to the rational fractions of n we have dealt with, we cannot conclude that 71 (q) = 3/2 
(the 3/2 comes from the Catalan number asymptotics) identically for large q. We expect 
this to hold only in the q 00 limit, whereas for finite q, 71(g) is some function of q. 
We can now write 



mM ^ coiq) Roiq)'' + c±{q) 



RijqT 



+ 



(5.4) 



where ± is chosen according to the parity of n and we have indexed by the functions 



R{q) and c{q) corresponding to the polynomial contributions to the M, 



(n-fc) 



with the 



asymptotic expansions (|4.12|) . When q is large, the first term dominates as Ro{q) ~ q, 
whereas Ri{q) ~ 2y^. This justifies a posteriori the identification R{q) = Ro{q) in the 
large q regime of sect. 4. The properties derived in sect. 4. for the large q regime will 
however break down if Ri{q) > RQ{q), i.e. when the corrections become dominant. We 
expect such a crossing of phases to take place at the transition point q = qc, ai which 
Ri{qc) = Ro{qc)- In this scenario, R{q) = Ri{q) for all q < qc, and the exponent 7(g) 
jumps abruptly from its large q > qc value 70 (q) = 0, to a non-zero value 7i(q'c) 7^ 0. This 
mechanism explains the jump in the value of 7(5) at the vicinity of q = qc, as observed in 
Fig.ll. In this figure, we have represented in dashed line the purported value 7(13') =0 for 
q > 2, which must be substituted to the (bad) large n estimate in this regime, whereas we 
still rely on the (good) q < 2 estimate. 
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To reconcile this scenario with the picture described in sect. 3, in which for q < Qc 
the winding of semi-meanders becomes irrelevant (i.e. R{q) = -R(<?)), we should have 
Ri{q) — R{q) below the transition point q^ where the semi-meanders enter their meander- 
like phase. Since the meanders themselves do not display any transition at q = qc, it is 
tempting to infer that [q) = R{q) for all values q. As we will see now, this is corroborated 
by the large q expansion ( |5.3| ) for Ri{q), which turns out to coincide with that of R{q) 
for meanders (|2.10|) . The latter is easily carried out for the meander numbers, whose 
structure, given in appendix E, is very similar to that of the semi-meander corrections. 



More precisely, the meander numbers Mn 



-k) 



take the general form 

Cn rk{n) 



(5.5) 



where r/j(x) is a rational fraction of x, with total degree k. With the values listed in 
appendix E, we find the following large q expansion for the functions R{q) and c{q) of 
( p.lOp , up to order 6 in 1/q 




(5.6) 

The exponent (y.{q) is also found to tend to the Catalan value 3/2 when q tends to infinity, 
but appears to be a non-constant function of q for all finite q. Remarkably, the first terms 
of the large q series expansions of R{q) and those of the correction to semi-meanders Ri{q) 
coincide! 

In conclusion, all our results conspire to suggest that the semi-meanders undergo at 
q = qc a, transition from a "meander" -like regime q < q^ governed by the meander partition 
function per bridge R{q) = Ri{q), to another regime q > qc, governed by Ro{q). The order 
parameter for this transition is clearly 



lim —{w)^ 

n— »oo n 



1%) 



for q > qc 
for q < qc 



(5.7) 



which vanishes for q < qc (irrelevant winding, i.e. u{q) < 1) and is nonzero for q > qc 
(relevant winding, i.e. u{q) = 1). With the order parameter (|5.7|) , the transition is found 
to be continuous, as the leading coefficient X{q) ( |4.18| ) vanishes at q = qc- The smooth 
character of the transition is also visible from the fact that R{q) and R{q) approach each 
other tangentially at q = qc, and that the coefficient Co{q) = R{q)a{q)/q of the large-g 
dominant contribution to rhniq) vanishes at q = qc (see Fig. 15). 
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6. Small q behavior of the semi-meander polynomial 



The very existence of asymptotics of the form ( p.6| ) for the semi-meander numbers, 
with a smooth enough function R{q) has highly non-trivial consequences on the numbers 
. We have already seen how the numbers , for large n and finite /c, are linked to 



each other ( [4 .81 ) in order for (|2.6|) to hold for large q. Let us now examine its consequences 



on the small q and large n behavior of mn{q). Let us expand rriniq) around g = up to 
order k, and take the large n asymptotics of each term in the expansion: 



m, 



n 

,2^ ' 



(6.1) 



k 

,{q) = Y.Mi^\^ +0{q^+^) 

3 = 1 

= 4-i?(o)-fi + ,£M + o(,^ 

n^ii) ^ ' \ R{0) ^ 

where we have only retained the leading n asymptotics in each q^ term, and used the 
q ^ limits ( |277| ) (actually, we have assumed that q ~ ^/n). Comparing with the 
asymptotics (pTTI), we finally get 

valid for large n and finite k. This is actually very similar to the behavior of the Mn" '^^ ~ 
n^/k\ for large n and finite /c, and may be deduced from the main recursion relation 

^ (k) 

for semi-meanders as well. Indeed, the number MA of semi-meanders of order n with 
k connected components is obtained from that of order 1 (root) by {k — 1) applications 
of the process (II) (see Fig. 8), and (n — k) applications of the process (I), whereas the 
Mn^'* connected semi-meanders are obtained through the process (I) only. Apart from the 
relative combinatorial factor (^Z}) ~ n^~^ /{k — 1)! accounting for the {k — 1) choices of 
process (II) among the total of (n — 1) steps, we must consider that whenever a step (II) 
is chosen instead of a step (I), some freedom in the overall choice is lost. Eg. ( |6.2| ) tells 
us that this corresponds to an average factor of -R'(0)/^(0) per step (II) taken instead 
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of a step (I). We checked (|6.2|) numerically, by performing a large n extrapolation of the 
appropriate ratio for a few values of k. We find a very good agreement with the estimate 



7. Conclusion 

In this paper, we have analyzed the meander problem in the language of critical phe- 
nomena, by analogy with Self-Avoiding Walks. In particular, we have displayed various 
scaling behaviors, involving both scaling exponents and scaling functions. We have pre- 
sented strong evidence for the existence of a phase transition for semi-meanders weighed 
by a factor q per connected component (road). In a large-g regime {q > qc), the winding is 
found to be relevant, with a winding exponent u{q) = 1, while the configuration exponent 
7(13') = 0. Throughout this phase, a finite fraction a{q)/X{q) of the winding is due to 
circles, i.e. circular roads with only one bridge, winding around the source of the river. In 
this regime, the partition function per bridge for semi-meanders R{q) is strictly larger than 
that of meanders R{q)- The particular form of its large q series expansion in 1/q ( [4.12| ) 
with slowly alternating integer coefficients, which furthermore grow very slowly with the 
order, and its purported re-summation (D.14), suggest a possible re-expressionion in terms 
of modular forms of yet to be found. In a low-q regime q < qc, R{q) and R{q) coincide, 
in agreement with an irrelevant winding i/[q) < 1. The exponent 7(g) is no longer 0, but 
a strictly positive function of q. 

We have estimated the value of the transition point qc — 2 with an excellent precision, 
and we conjecture that qc = 2 exactly. This special value of q has already been singled 
out in the algebraic study of the meander problem, in connection with the Temperley- 
Lieb algebra [0. There, we have been able to re-express the meander and semi-meander 
partition functions as that of some Restricted Solid-On-Solid model, whose Boltzmann 
weights are positive precisely iff g > 2, indicating very different behaviors for q < 2 and 
q>2. 

In the large-g phase, in addition to the exact values ^{q) = 1 and 7(13') = 0, we can use 
the asymptotic expansion (|4.12| ) to get R{q) with a very good precision. The somewhat 



sub-leading meander quantities R{q) and a{q) are more difficult to evaluate in this regime. 

A number of questions remain unsettled. There still remains to find the varying 
exponents 7((j) and ^{q) in the q < 2 regime, as well as the precise value of R{q) = R{q)- 



31 



Although we improved our numerical estimates, we are limited to conjectures. For q = 0, 
we confirm a previous conjecture that 7 = 2, and that a = 7/2. We also conclude 
from the numerical analysis that u{0) ~ 0.52(1) is definitely not equal to the trivial random- 
walk exponent 1/2. For q = 2, we have an amusing guess for R{2) = R{2) = 7vV2 = 
r(l/4)r(3/4) = 4.442... inspired from possible infinite product formulas for R{q). 
Acknowledgments 
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Appendix A. Algorithm for the enumeration of semi-meanders 

The following is a simple computer algorithm directly inspired by the recursion relation 
of Sect. 3.1. 

There, it has been shown that one can construct a tree of all the semi-meanders 
generated recursively with processes (I) and (II), as displayed in Fig. 8. Each node at depth 
n represents a semi-meander of order n. To have a finite cost of computation, the order is 
limited to nmax and the nodes of depth nmax appear to be leaves. The algorithm we used 
consists in visiting all the leaves from left to right, following the branches, as a (clever) 
squirrel. 

To do that, the rules are : 

(a) The squirrel starts at the root (upper node). 

(b) When the squirrel is on a intermediate node (not a leaf), he follows the leftmost branch 
which it has not yet visited and the depth increases by 1 . If all the branches of a given 
node have been visited, the squirrel goes back up one level and the depth decreases 
by 1. 

(c) When the squirrel is on a leaf (depth nmax), it goes back up one level and the depth 
becomes [nmax — 1). 

The reader can convince himself that the above rules describe a systematic and com- 
plete visit of the tree. Of course, when the squirrel is on a node, it can measure a lot 
of interesting quantities like the number of connected components, the winding number... 
These measures are added up and analyzed at the end of the enumeration. 

From a Fortran point of view, many representations of (semi-)meanders are possible. 
In the open-river formulation, each semi-meander is made of a lower rainbow arch config- 
uration (which we need not code) and an upper arch configuration of order n. For conve- 
nience, we label the 2n bridges of river from i = —n + 1 to z = n, and the system of arches 
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is described by a sequence of integers {A{i);i = —n+1, ...,n}, where A{i) e {— n + 1, n} 
is the label of the bridge connected with the bridge i. The following Fortran program enu- 
merates the connected semi-meanders. For simplicity, only the process (I) is coded and the 
number of connected components is always k = 1. The arch to be broken {j,A{j)) begins 
at the bridge j and ends at the bridge A{j). The process (I) splits this arch into two arches 
{—n,j) and {A{j),n+ 1). When the squirrel climbs back up one level, the two extremal 
arches {—n+ 1, A{—n + l)) and (A(n), n) are re-sealed to give one arch {A{—n+ 1), A{n)). 
At this stage, we know that the next arch to break starts at bridge j — A{n) + 1. This 
ensures the completeness of the algorithm. 



PARAMETER (nmax = 14) 
INTEGER A(-nmax+l:nmax) 
INTEGER Sm(nmax) 
INTEGER n 
INTEGER j 

DATA n, Sm /O, nmax*0/ 
A(0J = 1 

A(l) = 
n = n + 1 
Sm(n) = Sm(n) -|- 1 
j = -n + 1 

IF((n.EQ.nmax).OR.(j.EQ.n+l)) GOTO 3 

ArA(j)) = n+1 
A(n+1) = A(j) 
A(j) = -n 
A(-n) = j 

2(]5S+1)) = A(n) 
A(A(n)) = A(-n+l) 
j = A(n)+1 

fF"(n .GT. 1) GOTO 1 

PRINT '(i3, il5)', (n, Sm(n), n = 1, nmax) 

END 



maximal order 
arch representation 
semi-meander counter 
current depth (or order) 
next branch to visit 
n and Sm initialized to 
single- arch semi-meander 

! a new node is visited 

leftmost (exterior) arch 

up or down ? 

go down with process (I) 



! going up 

! next arch to break 



It is possible to use the left-right symmetry to divide the work by two. It is also 
possible to adapt the program for a parallel computer. For that, an intermediate size 
(nl = 11, for instance) is chosen. A first (little) run is made with nmax = nl, which gives 
M(nl) leaves. In a second (big and parallelized) run, each of these leaves is now taken 
as the root of a (sub-)tree and treated independently of the others. At the end, all the 
results of the sub-trees are collected. The calculations of this article have be done on the 
parallel Cray-T3D (128 processors) of the CEA-Grenoble, with approximately 7500 hours 
X processors. 
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We have computed Mn (the (pure) semi-meander number of order n) up to n = 29, 

(semi-meander number of order n with k connected components) up to n = 27 

and the other quantities up to n = 24. The reader can obtain an electronic copy of the 

numerical data upon request to the authors. We content ourselves with giving, on Table I 

— — 
below, the Mn and our last row (n = 27) of the M^'. 



n 


Mr, 


n 




k 


Z ( 


k 




1 


1 


16 


1053874 


1 


369192702554 


16 


2376167414 


2 


1 


17 


3328188 


2 


2266436498400 


17 


628492938 


3 


2 


18 


10274466 


3 


6454265995454 


18 


153966062 


4 


4 


19 


32786630 


4 


11409453277272 


19 


34735627 


5 


10 


20 


102511418 


5 


14161346139866 


20 


7159268 


6 


24 


21 


329903058 


6 


13266154255196 


21 


1333214 


7 


66 


22 


1042277722 


7 


9870806627980 


22 


220892 


8 


174 


23 


3377919260 


8 


6074897248976 


23 


31851 


9 


504 


24 


10765024432 


9 


3199508682588 


24 


3866 


10 


1406 


25 


35095839848 


10 


1483533803900 


25 


374 


11 


4210 


26 


112670468128 


11 


619231827340 


26 


26 


12 


12198 


27 


369192702554 


12 


236416286832 


27 


1 


13 


37378 


28 


1192724674590 


13 


83407238044 






14 


111278 


29 


3925446804750 


14 


27346198448 






15 


346846 






15 


8352021621 







Table I: The numbers MA of semi- meanders of order n with k connected components, 
obtained by exact enumeration on the computer: on the left, the one-component semi- 
meander numbers {k = 1) are given for n < 29; on the right, n is fixed to 27 and 1 < k < n. 

~ (k) 

The MA^ for ?i< 27 can be obtained by request to the authors. 



Appendix B. Correction terms for semi-meanders with large number of con- 
nected components 

In this appendix, we show that the corrections to the polynomial expression (^]^) for 
the numbers Mn" '^^ occur only for n < 2/c — 2, and derive the first correction fx^"^ of ( [4.6| ) 
for n = 2k — 2. The result reads 

Z^!,"^ = Ck-i (B.l) 
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where the are the Catalan numbers ( p.4| ). 

As explained in Sect. 4, the polynomial part of Mn is generically obtained as a 
sum of combinatorial factors, counting all the possible occurrences of perturbations of 
the leading semi-meander (order n, n components) which have the same order n, but 
have only n — k connected components. A perturbation, by definition, is made of a core, 
which does not contain any circle, supplemented by circles. The core is generically made 
of p irreducible semi-meanders with a total order no, and a total of (no — k) connected 
components. By irreducible, we mean that no circle can separate the semi-meander in two 
disconnected pieces. This core is completed by (n — no) circles, to form a semi-meander 
of order n. Enumerating all perturbations of the leading meander with n — k components 
amounts to enumerating all the ways of completing cores by circles. There are exactly 



n — uq + p 
p 



(n - no +p){n - uq + p - l)...(n - no + 1) 



fB.2) 



possible decorations of the above core by circles. The polynomial form (|B.2|) of the com- 
binatorial factor is however only valid for n > uq — p. When n < no — p — 1, we have to 
add a correction to the polynomial form to get a vanishing result, equal to (— 1)^+^ for 
the largest n = uq — p — 1. The largest n at which such corrections occur is obtained by 
maximizing uq — p. As the core does not contain any circle, each of its tt-q — /c connected 
components is at least of order 2, hence 



no > 2(no — k) 



no < 2k 



(B.3) 



This inequality is saturated for a core made of k kidneys. Minimizing p consists in taking 
{k — 1) arbitrary kidneys included in one, for which p = 1. We therefore find that the first 
correction to the polynomial behavior (|4.6D occurs at n = 2/c — 2. 





Fig. 16: Arbitrarily disjoint or included kidneys are equivalent to symmetric 
meanders upon folding back the river, as indicated by the arrows. The latter 
are in one-to-one correspondence with arch configurations of the same order. 
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There are exactly Cfc-i possible choices of these (/c — 1) kidneys, in one-to-one corre- 
spondence with arch configurations of order (/c — 1) as illustrated in Fig. 16, upon a folding 

(k) 

procedure. Each corresponding core contributes 1 to the correction /Xq , which completes 
the proof of (Q. 

More generally, we expect corrections to the polynomial part of MA for all k < 



n < 2/c — 2, hence the form (|4.6|) . In appendix C, the structure of the first 9 successive 
corrections is found, together with the polynomial part of MA , up to /c = 18. 



Appendix C. Fine structure of the semi-meander numbers 



As shown in appendix B, the polynomial part in ( [4.6|) of the semi-meander numbers 
Pk{n) is equal to the semi-meander number Mn" '^^ for n > 2k — 1. In addition, the 



polynomials pk are subject to the quasi- recursion relation (|4.9| ), which leaves only two 
coefficients of pk to be fixed, once the pi, I < k — 1 are known (we have po{n) = 1). It is 
therefore straightforward to derive the first 13 polynomials pk from our numerical data, 
as only the two values Pki^k — 1) = m!^_^ and Pki^k) = are required, and the 

corresponding numbers are known up to /c = 13. 

However, knowing the pkS up to /c = 13, we get a list of all the corrections 

ixf = M^lztf - Pk{2k - 2 - j) (C.l) 

to these polynomials needed to get the correct values of up to /c = 13. The 

resulting table of corrections displays a remarkable structure (very close to that of meander 
numbers, studied in appendix E) , which reads as follows, according to the parity of j 



= ^^7r7^X^^n4.-i(n) withn = 2p 



(p+i+j) _ ^ 2^' (2p+l) ' 
fit+i"'^ = ^p3|^^^^n4,(n) withn = 2p-l 



valid for all j > 0, where Cp are the Catalan numbers ( |2.4| ), and the nfc(n) are monic 
polynomials of degree k. 
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Proceeding in parallel (determining the pk and the corrections simultaneously), we 
can proceed further up to A; = 18, with the result 

Po{x) = 1 
Pi{x) = x-1 
P2ix) = {x'^ + x-8)/2 
psix) = {x^ + - 31a: - 24) /3! 
P4{x) = (x^ + 14x^-49x2 -254a;)/4! 
P5{x) = (a;^ + 25a;^ - 15a;3 - 1105a;2 - 1066a; + 1680)/5! 
peix) = {x^ + 39x^ + Ubx"^ - 2895a;3 - 10226a;^ + 8616a; + 31680)/6! 
P7{x) = {x^ + 56x^ + 532x^ - SllOa;^ - 50141a;^ - 201463;^ + 377208a; 
+ 282240)/7! 

ps{x) = (a;^ + 76a;^ + 1274a;^ - 5264a;^ - 165991a;^ - 422156a;3 + 19791 16a;^ 

+ 6031824a; + 1128960)/8! 
P9{x) = (a;^ + 99a;^ + 2526a;^ + 2646a;^- 41351 la;^- 2570589a;^ + 4826744a;^ 

+ 55185444a;2 + 54007920a; - 29756160)/9! 
Pio{x) = (a;^° + 125a;^ + 4470a;^ + 30090a;^ - 803607a;^ - 10282755a;^ 

- 6206320x'^ + 302065660x^ + 838000656x2 _ 17932032OX - 914457600)/10! 
pii{x) = (a;^^ + 154a;^^ + 7315a;^ + 96360a;^ - 1170477a;^ - 31531038a;^ 

- 1 167481 15a;^ + 1085347340a;^ + 7183991276a;3 + 4813856784a;2 

- 20917209600a;- 15487718400)/11! 

Pi2{x) = (x^2 + 186a:^^ + 11297a:^° + 231330x^ - 921657a:^ - 78859242a;^ 

- 632084629a;^ + 2301195270a;^ + 41279402956a;^ + 93554770056a;2 

- 181951879968a;2 - 528315782400a; - 281652940800)/12! 

p^^{x) = (x^^ + 221a;^2 ^ 16679a;^^ + 478621a;^° + 1380093a;^ - 164767317a;^ 

- 2382680443x^ - 496896257x*^ + 172979664286,T^ + 862378655996a;'^ 

- 433580004936a;^ - 7916932037664a;2 - 9843473134080a; 
-4807260057600)/13! 
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Pi4{x) = (x^^ + 259a;^^ + 23751a;^2 + 899171x^^ + 8661653x^° 

- 285047763a;^ - 7202313547a;^ - 30034731727a;^ + 538444152246^^ 
+ 5357996127484a;^ + 6470339335096a;^ - 6769757951 1744a;^ 

- 191470854038400a;^ - 156970678394880a; - 54050540544000)/14! 
Pi5{x) = {x^^ + 300a:^^ + 32830^^^ + 1575210^^^ + 26278252x^^ 

- 378408030a;^° - 18419182210a;^ - 164903537370a;^ + 1145711810243a;^ 
+ 25071998561610a;^ + 89293562501780a;^ - 342807873 156840a;^ 

- 2342105695034496a;3 - 3394586375786880a;2 - 1602031548902400a; 
+ 1953665505792000)/15! 

pie{x) = (a;^^ + 344a;^^ + 44260a;^^ + 2614640a;^^ + 63321622x^2 

- 244465312a;^^ - 40783574260a;^° - 621949980080a;9 + 743774155553a;^ 
+ 92574072382792a;^ + 650521234967240a;^ - 658128229828160a;^ 

- 19409328228712176x^ - 53224642530877824x^ - 35515613029674240x2 

+ 48894886046361600a; + 121101107871744000)/16! 
pnix) = (a;^^ + 391a;^^ + 58412a;^5 + 4155820a;^^ + 134381702a;^3 

+ 605876362a;^2 _ 78440970036a;^^ - 1905575960860a;^° 

- 7182404321807x^ + 271696014206903x^ + 3438145428617536x^ 

+ 5440767496706360a;^ - 113308202712264496a;^ - 594837570662080656a;^ 

- 742844861844630912a;3 + 780842535860839680a;2 + 4213612687918233600a; 
+ 2696110704967680000)/17! 

pisix) = (s^^ + 441a;^^ + 75684x^^ + 6372756a;^^ + 261854502x^^ 
+ 3203791542a;^^ - 128136377252a;^2 _ 5023361538468a;^^ 

- 45509495478447a;^° + 602361827803593a;^ + 14435051961445752a;^ 

+ 67015410007677768a;^ - 451434378729887216a;^ - 4879430147272561776a;^ 

- 11678841781394909184a;^ + 6041803064009266944a;^ 
+ 78098750156766044160a;2 + 148841653993843507200a; 

- 4417637856952320000)/18! 

(C.3) 
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The first II's read, 



^ / X o 33 2 49 

1397 13229 
n7(a;) = + 44a;^ + -^x^ + 3371a;^ + -^x^ + 11735a;^ + 6336a; + 1440 

nii(x) = x^i + l^x^V 2814^9 + ^^^x8 + 405471x^ + 1586409a;« 

+ 4522625a;^ + x"^ + 13577913x3 + 143219973;^ + 3974616a; 

4 

U,,{x) = a;i^ + 129x^4 + 7462x^3 + 248370x12 + ^^il^xii + iH^^a;i° 

2 2 

o.^oo^^o. Q -.o^^-...^^. 8 12479885017 7 44260908189 . 

+ 345882086a;^ + 1392155925a;^ + x^ + x^ 

2 2 

+ 30440330807a;^ + 71900294130x^ + 17319985020^2 

- 28987682544a;2 - 36059230080a; - 10059033600 
Ho (a;) = 1 

47 37 

U4{x) = x^ + —x^ + Slx"^ + —X + 20 (C.4) 

o o 

^ , , o 124 7 1283 18479 . 8543 4 72077 o 
Usix) = x^ + —x^ + ^x"" + —^x^ + ^x^ + —^x^ 

+ 26106a;2 - 4994a; - 840 

Ui2ix) = a;i2 + 77a;ii + ^^a;i° + ^^^a;^ + 341193a;^ + 1023192a;^ 

5147969 35955101 _ 

2 2 
+ 18500028a;2 - 21290040a; - 6350400 

^ , . 16 368 15 14 664235 13 8557681 12 135562460 n 

Hiefa;) = x^^ -\ a;^^ + 6845a;i^ H a;" H x^^ -\ a;" 

''^ ^ 3 3 2 3 

452366205 ^ 1940546440 „ 11216097197 o 89505296956 7 
H ^ a;^" H ^ x^ H ^ a;^ H ^ x^ 

^ 9305821085^, ^ 27505324345^, ^ 164921583110.^ - 422241984828.^ 

2 3 

- 328334271840a;2 - 252250683840a; - 40475635200 



Appendix D. Semi- meanders and winding 

In this appendix, we study the numbers Mn''\w) of semi-meanders of order n, with 
k connected components and winding w. The winding w has the same parity as the order 
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n, hence we consider the following generating function 

n-l min(fc,[n/2]) 

mn{q,t) = Y.q^-'^ J] Mi^-'^\n - 2j) (D.l) 

fc=0 j=0 

interpreted as the partition function for semi-meanders weighed by q per connected compo- 
nent and 1/Vi per winding unit (up to a global normalization factor). In the large q limit 
we now repeat the analysis performed in sect. 4.1, by Taylor-expanding Log{rhn{qTt) /q^) 
order by order in 1/q, in the same way as we did previously for Log{rhn{q) / q^) ■ This also 
relies on the identification of the numbers Mn"^ '^^ {n — 2j) as polynomials of n, with special 
re-exponentiation properties, as discussed now. 

For k = 0, the corresponding leading semi- meander has winding w = n, hence 

Mi^\n-2j) = 5,-o (D.2) 

For large n and finite k > 1, the winding of a semi- meander of order n with n — k connected 
components may only take the values (n — 2A;), (n — 2/c + 2), ...,(n — 2). Indeed, if C denotes 
the number of circles of such a semi-meander with winding {n — 2j), there are [n — k — C) 
connected components which are not circles, hence which occupy two or more bridges. 
Therefore we have the following lower bound on the order of the semi-meander 

n > C + 2{n-k-C) ^ C > n - 2k (D.3) 

and therefore, as each circle contributes 1 to the total winding, 

w > C>n-2k (D.4) 



When /c = 1, all the one-kidney perturbations (|4.3| ) of the leading semi-meander have 
winding w = {n — 2), hence 

M^"-^)(n-2) = n-l (D.5) 

When /c = 2, let us reexamine the various semi-meanders obtained in Fig. 14: the 
perturbations of Figs. 14 (c)i,2 have winding (n — 4), whereas those of Figs. 14 (0)3,4 have 
winding (n — 2). Hence we have 

M(-2)(n-4) = 

" ^ ^ 2 ""'^ (D.6) 

M^'^-2)(n-2) = 2(n-2) 
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The correction is ad hoc to yield a zero answer when n = 2. 

With a httle more patience, the enumeration of the semi-meanders of order n with 
(n — 3) connected components and fixed winding yields 

*("-»)(„ -6) = "'"-^^("-^' +2^3 + 2^,.,, 

M^"-3)(n-4) = 2(7i2 -4?! + 1) +45n,3 
= 2(n-3) 

In general, the form a decomposition of the ^\ in the sense that 

"^M^^-^\w) = M^"-'^) (D.8) 



In a way similar to the ^\ we expect the numbers — 2j) to be, for large 

enough n, some polynomials of n, whose coefficients depend only on j and k. For small n, 
some corrections have to be added to recover the actual numbers from their polynomial 
part. More precisely, one can show that 

i-2 

m^-^\n - 2j) = pf\n) + 5.,.+,-^-2 ryL^'^') (D.9) 

m=0 

for 1 < J < /c < n, where {x) is a polynomial of degree j of x {j>"q\x) = 1), whose 
coefficients depend on j and /c, and the rjrn are non-negative integer corrections. In 
particular, according to ( p.8| ), we must have 

fe 

Y,pf{x) = Pk{x) (D.IO) 
j=i 

We have computed the polynomials pj^"^ (x) for < j < k < 14, by using exact enumeration 
data on the Mn"^ ^""(n — 2j) for 0<j</c<n< 24, and the re-exponentiation trick 
described in sect. 4.1. These lead to the large n asymptotics of the partition function ( p.l| ) 
in the large q regime 

mn{q,t) ~ c{q,t)—:^^ (D.ll) 
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with 7(0', t) = as before, and with the following large q series expansions 



Riq,t) = q + t+- + - + -+ ^ ' +^ 

q q-^ q* q^ 

2t{2,t^ - 8t + 1) 2t{W^ - St^ + - I2t + 1) 
fi 7 

2t{10t^ - 26t^ + 28t^ - 18t + 1) 2t{6f^ + 12t^ - 61t^ + 64t^ - 24t + 1) 

^ ^8 ^ ^9 

2t{24t^ - 52t^ + 71t^ - 137^3 + 131t^ - 32t + 1) 
+ ^10 

2t{101t^ - 26Qt^ + 308t^ - 324t^ + 237t^ - 40t + 1) 
+ q^' 

2t{90f + 50t^ - 610t^ + 894t^ - 726t^ + 400t^ - 50t + 1) 

+ ^12 

^ 2t(173t8 - 243t^ + 409^6 - 1507t^ + 2237t^ - 1564^^ + 631^2 - 60t + 1) 



qi3 



(D.12) 



, , t At 2t(t-S) 8t(t-l) 2t(7t'^ - m + 5) 

c{q,t) = 1 2+ 3 + 4 5 

q q^ q-^ q^ q° 

2t{2t^ - 13*2 + 39t - 6) 2t{18t^ - 62*^ + 79t - 7) 

^ ^6 + ^7 

2t(45t^ - 122^3 + 183^2 - 128t + 8) 

2t{9t^ - 138t^ + 406^3 - 514^2 + 205t - 9) 

+ q' 

2t(69t^ - 475^4 + 1138^3 - 1143^^ + 295t - 10) 
+ ^10 

2t{353t^ - 1183^5 + 2104t^ - 3012*^ + 2342t2 - 421t + 11) 

2t(56t' - 1553t*^ + met'' - 7692^4 + 7668f'^ - ISllf^ + 56U - 12) 
+ qi2 

_ 2t(230t^ + 3888^6 - 14710t^ + 22124t^ - 17898*^ + 7490t2 - 751t + 13) 

^13 

2f 

- -^(3123t8 - 10053t^ + 22655t^ - 45694t5 + 58185*^ - 39176*^ 
+ 12215^2 - 959t + 14) + 0(4f) 

(D.13) 
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The previous results ( [4.12| ) are recovered up to order 14 in by taking t = 1. Let us look 
more closely at the expressions ( p.l2|) and ( p.l3| ) above: we may Taylor-expand them as 
functions of t for small t. Remarkably, due to the structure of the coefficients of pf\x), 



which are themselves polynomials of k for fixed j (e.g. we have p^i\x) = 2{x — k) for all 
k > 2 while p\^\x) = x — 1), we have been able to re-sum the large q series coefficients of 
this expansion up to order 3 in t, in the following form 

R{q^t) =q + t- -At^ 



q-l (Q-l)2(g2_i) 
, .3 - + + llg^ +8q^ + q-l . 

Clg,tJ - g(g-l)2 1)2(^2 _ 1)2 

- %4(, _i)(g2! 1)3(^3 _ 1)2 (V^ + 6,^^ + 28,^'^ + 8V + 193g« 



+ 228g^ + 147^^ + 35^^ - 14^^ - 5q^ + Gq"^ + q - 1) + 0{t 
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This structure is very reminiscent of the large Q Q-states Potts model free energy ||Tl 
and suggests a possible expression in terms of modular forms of q. 

The average winding in large q semi-meanders is obtained by the formula 

{w)niq) = n-2t^Logmn{q,t)\^^^ (D.15) 
leading to the expression ( [4.17| ) through the identifications 

X{q) = l-2t-^Log^(Q,t)|^^^ 
M = -'^t— Log c{q,t)\ 



\t=i 



Appendix E. Fine structure of the meander numbers 

By simple inspection of the meander numbers (which we read from the semi-meander 
numbers with zero winding), we have found the following structure for the numbers '^^ 

M^-'^) = nin) (E.l) 

where is a Catalan number ( ^.4] ), rk{x) is a rational fraction of x with total degree k. 
More precisely, we have, ro{n) = 1 and for k > 1 

2"n!(n + l)! 
''^"^ = fe!(n-fc-l)!(n + 2fc)! ^^^-^^"^ 
2'' n{n-l){n-2)...{n-k) 
¥ (n + 2)(n + 3)...(n + 2k) ^^fc-sl^J 
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where ^2k-2{x) are monic polynomials of x of degree {2k — 2). They read, for A; = 1, 2, 6 

(Pq{x) = 1 

(p2ix) = x'^ + 7x — 2 

(P4{x) = a;^ + 20a;^ + 107a;^ - 107a; + 15 

m^fx) = a;^ + 39x^ + 547x'* + 2565,T"'^ - 5474x2 + 2382x - 672 

(E.3) 

ips{x) = a;^ + 64a;^ + 1646a;^ + 20074a;^ + 83669a;^ - 323444a;^ 

+ 257134a;2 - 155604a; + 45360 
(fiio{x) = a;^° + 95a;^ + 3840a;^ + 83070a;^ + 940443a;^ + 3382455a;^ 

- 22294735a;^ + 27662860x^ - 26147139a;2 + 16354530a; - 4098600 
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